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Optical  Matrix  Inverter  for  Phased  Array  Radar 


I.  Introduction 

The  purpose  of  this  study  is  to  investigate  new  approaches  to  process  the  data 
collected  by  phased  array  antennas  to  achieve  interference  canceling.  This  data  processing 
need  to  be  performed  with  severe  time  constraints.  This  study  was  motivated  by  the  fact 
that  optics  is  capable  of  processing  data  in  parallel  with  very  high  speeds.  The  issue  we 
investigated  in  this  study  was  the  applicability  of  the  optical  algebraic  processors  to  such 
problems.  The  advantages  gained  by  using  optical  systems  over  their  electronic 
counterparts  are  investigated.  The  limitations  of  the  optical  system  and  their  effect  on  the 
interference  canceling  was  studied.  In  the  course  of  this  contract  we  investigated  the  issues 
just  listed  and  came  to  valid  and  encouraging  conclusions  which  are  outlined  in  the  final 
section  of  this  report. 

II.  Interference  Canceling 

Phased  array  antennas  are  used  both  in  radar  and  communication  systems  among 
many  other  applications.  These  systems  are  used  to  detect  a  desired  signal.  In  practical 
systems  the  desired  signal  will  not  be  the  only  signal  present.  Thermal  noise,  other  friendly 
and  unfriendly  interference  signals  usually  are  present  in  surrounding  environment.  It  is 
common  that  the  jamming  signals  be  stronger  than  the  desired  signal,  which  presents  a 
serious  problem  in  detecting  the  desired  signal.  It  was  suggested,  to  resolve  this  problem, 
that  if  we  can  change  the  antenna  pattern  in  such  a  way  that  we  introduce  a  null  along  the 
direction  of  the  jammer,  it  will  cancel  it  without  effecting  the  ability  of  detecting  the 
desired  signal  [1—5].  This  change  of  the  antenna  pattern  can  be  achieved  by  adapting  the 
weights  of  the  individual  sensors  to  the  environment.  This  adaptation  process  needs  to  be 
done  using  the  detected  signals.  This  is  the  basic  idea  behind  the  adaptive  phased  array 


antenna  systems.  In  the  following  section  we  show  how  this  adaptation  processes  can  be 
achieved. 

2.1  Adaptive  Phased  Array  Antennas 

In  adaptive  phased  arrays  the  incoming  signal  is  detected  by  an  array  of  sensors. 
The  detected  signal  is  a  combination  of  the  target  signal  plus  interference  and  noise  signals. 
The  system  is  adjusted  in  such  a  way  to  suppress  the  interference  signals  reception  without 
affecting  the  desired  signal. 

In  this  section  we  consider  the  two  general  cases  of  interference  canceling:  first  by 
assuming  that  the  interference  signal  direction  is  known;  secondly  by  assuming  no  a  prion 
information  is  known  about  the  interference  signal. 

A.  Interference  Signal  Direction  is  Known 

When  the  interference  signal  direction  is  known  the  weights  wj’s  of  the  array  can  be 
chosen  to  suppress  the  interference  signal.  Let  the  system  shown  in  Fig.  1(a)  be  used  to 
demonstrate  this  adaptation  technique.  The  output  signal  of  the  array  s(t)  is  given  by  1 
s(t)  =  P[(w1  +  w^)  sin  uQ t  +  ( w 2 + w ^ )  sin(a’Qt  -  0-^)} 

+  I[w^  sin(cuQt  —  0)  +  v/y  sin(wot  —  0  — 

+  Wg  sin(wot  +  0)  +  w^  sin(a;ot  +  0  —  ^)]  ,  (1) 

where 

P  =  the  pilot  signal, 

I  =  the  interference  signal,  and 
0  =  the  phase  shift 

0  =  sintA  (2) 

To  cancel  the  interference  signal  and  to  make  the  signal  s(t)  equal  to  the  pilot 
signal,  we  need  to  solve  the  following  system  of  linear  equations  for  the  weights  \w’s: 
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Figure  1  Basic  adaptive  array  system  with  (a)  signal  and  noise  directions  are  known, 
and  (b)  no  a  priori  information  is  assumed. 


(3) 


W1  +  w3  =  1 

w2  +  w4  =  0 

(w4  +  )  Costf  —  (w’2  —  w^)  sintf  =  0 

(v?2  +  )  Cost?  +  (w,  —  w^)  sin#  =  0 


The  size  of  this  system  of  linear  equations  depends  on  the  number  of  sensors  in  the 
array.  The  number  of  jammers  can  make  the  system  under  or  overdetermined,  which  are 
both  time  consuming  algebra  problems. 


B.  No  a  priori  Information  is  Known 

This  is  the  most  general  case  where  we  assume  no  information  about  jammers.  The 
system  used  in  this  case  is  shown  in  Fig.  1(b).  Each  of  the  n  sensors  receives  a  signal  x;(t) 
which  is  in  turn  multiplied  by  a  variable  weight  Wj.  The  output  signal  s(t)  is  compared 
with  the  desired  signal  d(t)  and  their  difference,  the  error  signal  e(t),  is  used  to  determine 
the  value  of  wj’s.  The  output  of  the  array  is 


s(t)  =  2,  Xj(t)  Wj 


i=l 


or 


where 


s(t)  =  2 


(4) 


(5) 


'v 

'  x^t) ' 

\t’  = 

w  ■ 

1 

ana  Z  = 

Xj(t) 

w 

i„(‘) . 

For  digital  sampled  data 


s(j)  =  wT  *(j)  , 


and 


<(j)  =  d(j)  -  s(j)  =  d(j)  -  i(j)  • 


(6) 


(7) 


(8) 


The  optimum  value  of  the  weights,  w-’s,  is  the  one  reduces  e(j)  to  zero  or  at  least  minimize 


For  N  samples  of  data  the  optimum  weights  satisfy  the  following  set  of  systems  of 
linear  equations: 

£(  1 )  =  d(  1 )  ' 

^  H  i  )  =  d(  i  )  (9) 

w^  i(N)  =  d(N)  .  _ 

The  N  sets  of  equations  have  n  unknowns,  and  usually  N  >>  n,  and  are  inconsistent 
and  over  specified.  The  optimization  problem  can  be  rewritten  as 


where 


vfr  ,  =  R  *  ?  , 
opt  xx  xd 


Rxx  =  E{“  }  > 


txd  =  E{5d}.  (12) 

The  matrix  is  called  the  covariance  matrix,  where  E{  • }  is  the  ensemble  average. 

Many  algorithms  are  introduced  [2]  to  solve  for  the  weights  in  Eq.  (10).  Some  of 
the  popular  algorithms  are  the  least  mean  square  (LMS),  and  the  direct  matrix  inversion 
(DMI). 

We’ll  briefly  mention  the  DMI  algorithm  since  it  leads  to  the  algorithm  introduced 
in  this  paper.  Eq.  (10)  cannot  be  determined  exactly  using  a  limited  number  of  samples  of 
the  input  data.  For  practical  consideration  a  small  number  of  samples  is  detected  to  be 
used  in  determining  w.  The  estimated  value  of  Eq.(10)  can  be  given  by 


W  =  R  1  ? 

xx  xd 


where 


Rxx  is  the  sample  covariance  matrix,  and  is  the  sample  cross— correlation  vector,  and 
are  given  by 


(14) 


K 

*xx  =  K  X  *(j)  iT(j) 

J=1 

and 

K 

?xd  =  K  X  ’  (15) 

'  j=l 

and  K  is  the  number  of  samples.  The  DM!  algorithm  determines  the  inverse  of  the  sample 
-♦ 

covariance  matrix  R  then  from  Eq.  (13)  evaluates 

III.  Optical  Algebraic  Processors 

Solving  systems  of  linear  equations  and  determining  the  eigenvalue  and  the 
eigenvectors  are  only  a  few  of  the  challenging  problems  faced  by  the  numerical 
computations.  The  problem  of  determining  the  weights  for  the  adaptive  phased  array  as 
given  by  Eq.  (13)  is  solving  a  system  of  linear  equations.  In  this  case  the  size  of  the 
matrices  involved  are  very  large.  Solving  a  system  of  linear  equations  for  large  matrices  is 
time  consuming  because  of  the  computation  complexity.  Digital  computers  revolutionized 
this  field,  because  of  the  fast  execution  of  number  crunching  operations.  But  still  solving  a 
problem  with  a  matrix  of  1000*1000  elements  takes  few  seconds,  which  by  our  standards,  is 
a  long  time. 

Optics  by  its  inherent  parallelism  and  speed  seems  to  present  a  natural  choice  for 
solving  this  class  of  problems.  Analog  optics  is  very  attractive  for  optical  information 
processing  and  computing.  As  shown  in  Fig.  2,  in  the  vector— matrix  multiplier  all  the 
elements  of  the  vector  are  processed  in  the  same  time.  At  almost  the  same  time  we  write  x 
we  do  read  b.  If  the  optical  path  length  between  the  input  and  output  planes  is  3cm,  the 
whole  operation  of  the  vector-matrix  operation  can  be  done  in  less  than  lOOpsec.  For 
X=1000,  the  number  of  operations  needed  to  perform  Ax— b  is  0 ( 1 0 6 ) .  Hence,  speed  of  the 
processor  is  0(10—  )  sec/operation.  This  illustrative  example  gives  a  sense  of  the  speed  of 


b 


J 


Figure  2  Incoherent  vector-natrlx  nultlplier 


the  analog  optics  in  performing  linear  algebra  operations.  Unfortunately  this  high  speed  of 
operations  is  combined  with  a  low  accuracy,  which  is  the  nature  of  the  all  analog  systems. 

Analog  optics  is  very  fast  but  inaccurate.  On  the  other  hand,  digital  electronics  is 
very  accurate  but  not  as  fast  as  analog  optics.  Utilizing  the  advantages  of  both  analog 
optics  and  digital  electronics  can  be  achieved  in  a  "compromise"  hybrid  system.  A  system 
that  slows  down  the  processor  speed  but  in  return  increases  the  accuracy  substantially. 

The  bimodal  optical  computer  (BOC)  introduced  by  Caulfield,  et.  al.  [6]  is  based  on 
this  idea  of  combining  the  speed  of  analog  optics  and  the  accuracy  of  digital  electronics. 

The  adaptation  problem  for  the  weights  introduced  in  Section  II,  can  be  rewritten 
in  the  following  form,  from  Eq.  (13) 


R  w  =  f  . 
xx  xd 

which  can  be  written  as 
AX  =  b  , 

where 


’16) 


:i7) 


A  =  R 


xx  ' 

X  =  , 

and  b  =  ?xd  • 

Eq.  (1G)  is  a  system  of  linear  equations  can  be  solved  using  the  bimodal  optical  computer. 

Consider  an  N*N  matrix  A,  and  N*  1  vectors  x  and  b.  Let  A  and  b  are  given,  and 
we  would  like  to  solve  the  system  of  equations  given  by  Eq.  (17)  for  the  vector  x.  This  can 
be  solved  by  analog  optics  techniques.  The  relaxation  method  introduced  by  Cheng  an 
Caul  Geld  [7]  can  be  used  to  solve  Eq.  (17)  for  x.  Consider  the  hybrid  system  shown  in  Fig. 

3.  Assume  an  initial  value  for  the  solution  x  and  write  it  using  the  LED's.  Then  the 

-*  * 

vector  x  is  multiplied  by  the  matrix  A.  The  resultant  vector  y  is  compared  with  b  by  a 

difference  amplifier.  This  difference  is  fed  back  to  correct  x.  This  process  of  multiplying 
the  new  value  of  x  with  A  and  comparing  y  to  b  continues  till  the  difference  between  y  and 


♦  -4 

b  becomes  zero  Then  the  value  of  x  will  converge  to  the  solution  of  Eq.  (17).  For  a 
positive  definite  matrix  A  always  a  convergence  to  the  solution  exists.  To  achieve  a 
nunnegative  definite  matrix,  we  can  multiply  Eq.  (17)  from  the  left  by  the  Hermitian  A11  of 
A.  The  new  A11  A  is  non— nonnegative  definite.  We  will  show  later  that  the  increase  in 
condition  number  this  causes  need  not  affect  convergence  and  that  we  can  .v.lve  the 
equations  even  if  A 11 A  is  singular. 

This  method  in  solving  a  system  of  linear  equations  is  very  rapid.  Its  speed  is 
limited  only  by  the  speed  of  the  electro— optics  devices  and  on  the  feedback  electronics, 
which  can  be  in  the  psec  range. 

Let  us  consider  now  the  accuracy  of  the  system.  In  writing  both  A  and  x  on  the 
optical  mask  (it  can  be  a  photographic  film  or  a  spatial  light  modulator)  and  the  LED 
array,  a  considerable  amount  of  error  will  exist  because  of  the  nature  of  these  analog 
devices.  Also  reading  the  vector  b  on  the  photodiode  array  cannot  be  done  exactly. 
Therefore,  the  system  in  Fig.  3  did  not  solve  the  system  in  Eq.  (17)  but  instead  the  system 
given  by 


A0x0  =  b0'  (18> 

where  the  subscript  zeros  indicate  inaccuracies  in  the  optics  and  electronics.  The  solution 
-*  -* 

Xq  of  Eq.  (18)  can  be  reGned  to  get  the  vector  x  using  the  following  algorithm: 


(a)  Solve  the  system  in  Eq.(18)  using  the  analog  optical  processor  for  Xq. 

(b)  Store  the  solutions  Xq  to  a  high  accuracy  with  the  digital  processor.  Use  a 
dedicated  digital  processor  to  calculate  the  residue 

r  =  b  —  Axq  =  A(x  —  Xq)  =  A  Ax  .  (19) 

(c)  Use  the  optical  analog  processor  to  solve  the  new  system  of  linear  equations 

Vo  =  sro-  (20) 

'4  -4 

where  y  =  s  Ax  and  s  is  a  "radix,"  or  scale  factor,  chosen  to  good  use  of  the  dynamic  range 
of  the  system. 


A 


Figure  j  The  bimodai  optical  computer  used  in  solving  a  system  of  linear  equations. 
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(d)  Use  the  digital  processor  to  refine  the  solution  Xq  for  x^: 

X1  =  XQ  +  Ax-  (21) 

If  the  refined  solution  x^  is  accurate  enough,  terminate  the  iterations.  Otherwise,  return  to 

(b),  (c),  and  (d)  for  a  more  refined  solution.  This  system  which  implements  the  algorithm 

outlined  above  is  shown  in  Fig.  3. 

The  convergence  and  speed  of  the  solution  for  the  system  of  linear  equations  is 
studied  and  reported  by  Abushagur  and  Caulfield  [8],  The  convergence  of  the  iterative 
solution  depends  on  two  main  factors.  First,  is  the  condition  number  of  the  matrix 

Aq,x(Aq).  The  smaller  condition  number  the  faster  it  will  converge.  Secondly,  on  the 

*♦  -♦ 

error  involved  in  reading  and  writing  of  A,  x  and  b  using  the  electro— optic  devices.  The 
higher  the  accuracy  in  representing  these  parameters  the  faster  the  convergence  will  occur. 

The  condition  number  of  the  matrix  is  a  critical  factor  in  the  convergence  of  the 
solution.  In  representing  the  matrix  A  by  an  optical  mask,  an  error  will  be  added  to  it. 
The  inaccuracies  in  representing  the  matrix  A  changes  the  values  of  the  matrix  elements 
a^’s.  These  variations  in  the  matrix  elements  change  the  condition  number  of  the  matrix. 
Let  us  represent  the  mask’s  matrix  Aq  in  the  following  form 

A0  =  A  +  E,  (22) 

where  E  is  an  error  matrix.  The  error  matrix  E  is  generated  using  Gaussian  statistics,  with 
standard  deviation,  a g. 

The  effect  of  the  error  matrix  E  on  the  condition  number  of  the  optical  mask’s 
matrix  is  demonstrated  in  Fig.  4.  In  Fig.  4(a),  a  matrix  A  with  condition  number  of  60  is 
considered.  The  coefficients  of  the  matrix  are  normalized  such  that  the  maximum  a-  is 
equal  to  the  unity.  The  condition  number  of  Aq  plotted  as  a  function  of  the  standard 
deviation  of  the  error  matrix,  The  condition  number  of  the  matrix  Aq  tends  to 

decrease  by  the  increase  of  <7g,  especially  for  large  In  Fig.  4(b)  a  matrix  A  with 
condition  number  300  (an  ill-conditioned),  the  condition  number  decreased  significantly 


throughout  the  range  considered  of  a  .  Thus,  if  A  is  an  ill-conditioned  matrix,  the  mask’s 
matrix  can  very  well  be  a  better  conditioned  one.  Of  course,  in  this  case  when  we  solve  the 
system  give  in  Eq.  (IS)  A(J  will  be  different  from  the  original  A.  Hence,  we  solve  a  better 
conditioned  system  for  the  approximate  solution  Xg,  and  then  we  refine  using  the  algorithm 
outlined  above. 

Now,  let  us  consider  the  effect  of  the  condition  number  on  the  convergence  of  the 
solution.  The  condition  number  is  the  determining  factor  in  the  accuracy  of  the  solution  of 
the  system  of  equations. 

Hence,  for  a  matrix  with  a  large  condition  number,  the  first  iteration  of  the  solution 
with  a  limited  accuracy  computer  will  be  highly  inaccurate.  This  leads  to  the  result  that 
the  larger  the  condition  number  the  larger  number  of  iterations  needed  for  the  convergence 
of  the  solution.  To  demonstrate  this  result,  we  ran  a  computer  simulation  of  our  bimodal 
optical  computer.  The  simulated  BOC  is  used  to  solve  a  system  of  linear  equation  with  a 
16  bit  resolution.  The  matrix  A,  was  generated  randomly  using  Gaussian  statistics.  An 
error  matrix  E,  with  an  error  of  1%  of  that  of  the  maximum  coefficient  of  the  matrix  A, 
and  then  added  to  A  to  generate  Aq  as  in  Eq.  (22).  An  error  of  1%  also  used  in  reading  Xq 
and  in  writing  bg.  In  each  case  we  computed  the  condition  numbers  of  the  matrix  and  its 
mask.  The  number  of  iterations  required  for  convergence  of  the  solution  to  the  specified 
accuracy  was  determined  for  each  case.  The  iterations  were  terminated  if  they  exceeded  25 
or  ||  r  (k+1^ll/l|r(k)||>  1,  which  is  the  condition  for  the  solution  divergence.  The  number  of 
iterations,  Nj,  required  for  convergence  of  the  solution  with  16  bit  accuracy  is  plotted  as  a 
function  of  the  condition  number  a(Aq)  *n  Fig.  5.  In  these  experiments  it  is  clear  that  the 
number  of  iteration  increases  with  the  increase  of  the  condition  number. 

The  condition  number,  as  shown  above,  is  one  of  the  determining  factors  for  the 
number  of  iteration  required  for  convergence  of  the  solution  of  the  system  of  equations.  It 
is  also  shown  in  Fig.  4  that  the  condition  number  of  the  optical  mask’s  matrix  decreases  by 
the  increase  of  the  standard  deviation  of  the  error  matrix  E. 
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5  The  nuaber  of  iterations,  N . ,  ia  plotted  as  a  function  of 


The  influence  of  the  standard  deviation  of  the  error  matrix,  E,  on  the  convergence  of 
the  solution  is  shown  in  Fig.  6.  'The  number  of  iteration,  Np  increases  with  the  increase  of 
jp.  This  decrease  in  the  convergence  rate  is  expected  because  for  large  a^’s  the  matrix  Aq 
is  quiet  different  from  A.  The  important  result  demonstrated  in  Fig.  6  is  that  even  with  an 
error  up  to  50%  in  writing  the  matrix  A  in  the  optical  mask,  convergence  is  still  achieved. 

The  result  is  very  important  in  realizing  this  algorithm  by  analog  optics.  In 
representing  the  matrix  A  by  an  optical  mask  always  an  error  will  exist.  An  error  of  1%  is 
quite  hard  to  achieve  in  this  representation  using  our  current  technology.  In  the  present 
state-of— the— art  technology  an  accuracy  of  2  to  3%  in  writing  the  matrix  A  is  within  our 
reach.  This  accuracy  does  not  sacrifice  the  convergence  of  the  solution. 

The  above  results  show  that  the  Bimodal  Optical  Computer  can  solve  a  system  of 
linear  equations  with  very  high  accuracy.  This  accuracy  can  be  achieved  using  I/O  devices 
that  have  limited  accuracy.  The  digital  computes  are  capable  of  achieving  high  accuracy 
solution  fur  all  the  cases  considered  above.  So,  what  is  the  real  advantage  of  introducing 
this  new  class  of  computers9  Speed,  is  what  we  are  after.  An  analysis  of  the  speed  of  the 
HOC  shows  that  for  it  to  be  more  faster  than  the  digital  computer  the  following  condition 
should  be  satisfied  [8; 


A  A,  >>  1, 
p  I 


where 


A  = 
P 


2[N°/6)  +  NT(1— K  -  NK] 


AI  -  TDpTAl  ’ 


K  =  I0/1D  , 


N  =  the  size  of  the  matrix  , 


(27) 


the  standard  deviation.  av.  of  the 


Ig  =  the  number  of  iterations  needed  for  the  convergence  of  the  solution  to  the  specified 
accuracy  using  the  BOC. 

Ip  =  the  number  of  iterations  needed  for  the  convergence  of  the  solution  using  the  digital 
computer. 

fpi  =  the  time  required  to  perform  one  digital  operation,  and 
T  ^  =  the  time  required  to  solve  AqXq  =  b,j  using  the  analog  processor. 

The  speed  advantage  depends  on  the  size  of  the  matrix,  X,  and  the  speed  of  the 
electronic  and  electro— optic  devices  used  in  the  BOC.  The  factor  A  is  plotted  in  Fig.  7  as 
a  function  of  X,  for  a  set  of  values  of  K.  It  is  quiet  clear  that  A  is  very  large  0(10°)  for 
moderately  large  values  of  X. 

The  values  of  Tp^  and  T  ^  can  be  compared  using  approximate  values  using 
current  data. 

T^  a  2  /isec  (27) 

Tpj  a  1  fisec  for  a  microcomputer,  (28) 

and 

TD1  a  1  nsec  for  a  CRAY2.  (29) 

The  factor  A  Aj  is  plotted  in  Fig.  8  as  a  function  of  N  using  the  data  given  by  Eqs. 
(27)  and  (29).  The  advantage  in  speed  is  very  large  and  the  condition  of  eq.  (23)  is 
satisfied  for  N>50. 

This  advantage  in  speed  of  the  BOC  over  the  existing  digital  computer  makes  it  a 
very  attractive  computing  machine,  and  shows  the  potential  of  this  class  of  hybrid  systems. 

IV.  Impelmenting  the  BOC 

The  BOC  was  built  in  our  laboratory  has  three  main  parts  as  shown  in  Fig.  3.  The 
optical  system,  the  electronic  circuit,  and  the  digital  processor.  The  optical  system  consists 
of  the  fully  parallel  matrix-vector  multiplier.  Light  from  the  LED’s  representing  x 
components  are  spread  vertically  by  planner  waveguides  onto  the  columns  of  the  matrix 
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mask.  The  transmitted  light  is  summed  row  wise  by  using  another  set  of  planner 
waveguides  and  detected  by  photodiodes  which  represent  the  output  vector  b. 

The  electronic  circuit  acts  as  a  feedback  loop  to  correct  for  the  input  light  of  the 
LED’s,  until  a  solution  is  reached.  The  solution  x  will  then  be  read  and  stored  by  the 
digital  processor.  Figure  8  shows  the  electronic  circuit  used  for  the  feedback  loop. 

The  A/D  and  D/A  conversion  from  and  to  the  electronic  circuit  are  performed  by 
the  digital  processor. 

In  this  section  we  present  the  experimental  results  for  solving  a  system  of  linear 

-♦  -♦  -♦ 

equations  Ax  =  b  using  the  BOC,  where  A,b,  and  x  are  ail  positive. 

The  Log  of  the  error  and  that  of  the  residue  are  plotted  versus  the  number  of 
iterations.  The  error  and  the  residue  are  defined  as, 

Error  =  ||(x-xk)||/||x||  (30) 

Residue=||  r^||  (31) 

Where  ||.||,  is  the  Enclidean  norm,  x  is  the  exact  solution,  x  is  the  k  iteration 
-»k  t  h 

solution,  and  r  is  the  k  iteration  residue. 

Since  we  are  dealing  only  with  positive  numbers  in  this  paper,  we  used  the  absolute 
value  of  r  to  solve  Eq.(2),  then  we  set: 

;(”+i)  =  +  (32) 

-4 

when  all  the  components  of  r  are  positive.  And 

+  =  {33) 

-»  ■* 

if  all  the  components  of  r  are  negative.  We  reject  the  iteration  when  the  components  of  r 

have  different  signs  and  take  the  previous  one.  By  rejecting  some  iterations  we  are  actually 

rejecting  some  corrections  which  will  slow  the  convergence  process. 

In  all  experiments  performed,  the  iteration  process  is  stopped  when  a  16  bit 

accuracy  is  reached.  Fig  9  shows  that  the  BOC  started  almost  30%  error  and  it  needed  6 


iterations  to  converge  to  1G  bit  accuracy.  In  Fig  10(a)  BOC  started  with  almost  1H)(% 
error,  and  the  number  of  iterations  needed  was  21.  Fig  10(b)  shows  the  Log  of  the  residue 
as  a  function  of  the  number  of  iterations.  The  fluctuations  depicted  by  Figs.  10(a)  and  (b) 
is  due  to  the  rejection  method  used  in  the  experiments. 

4.1  Effect  of  Calibration. 

The  analog  optical  system  error  is  a  major  factor  in  the  rate  of  convergence  of  the 
BOC.  If  that  error  is  reduced,  then  the  convergence  is  much  faster.  In  order  to  illustrate 
this,  the  same  problem  has  been  solved  twice  with  two  different  accuracies  of  the  optical 
system.  The  analog  optical  system’s  error  in  the  First  time  was  50%,  and  it  was  30c/c  in  the 
second  time.  Twenty  one  iterations  were  needed  by  BOC  to  converge  to  the  16  bit  accuracy 
for  the  first  case.  For  the  second  case  the  number  of  iterations  was  reduced  to  12.  These 
results  are  plotted  in  Fig  11. 

4.2  Reliability  of  the  System 

System  reliability  for  convergence  have  been  tested  and  verified  by  solving  the  same 
problem  several  times,  under  different  conditions.  Results  show  that  when  the  BOC  is 
used,  to  solve  a  problem  several  times,  the  convergence  rate  will  not  be  exactly  the  same 
for  all  the  cases.  However,  the  number  of  iterations  needed  by  the  BOC  to  converge  to  a 
certain  accuracy  is  almost  the  same.  Fig  12  shows  three  different  paths  of  convergence  for 
the  same  problem.  The  BOC  needed  13  iterations  in  the  first  run,  14  iterations  in  the 
second,  and  11  in  the  third. 

4.3  Convergence  of  the  Singular  Matrix 

Solving  a  system  of  linear  equations  with  a  singular  matrix  A  is  one  of  the  problems 
that  cannot  be  solved  using  conventional  digital  computer  techniques.  Singular  matrices 
have  a  condition  number  that  is  equal  to  infinity,  so  their  inverse  docs  not  exist,  also  they 
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Figiu(b;  The  Log(residue)  as  a  function  of 
the  number  of  iterations. 


Fig  10 ,  a )  The  Log(error')  as  a  function  of 
the  number  of  iterations.  The  BOC  started 
with  100%  error. 


Fig  1 1  The  Log(error)  as  a  function  of  the 
number  of  iterations  for  the  same  problem, 
but  with  two  different  accuracies  of  the 
optical  system. 


Fig !  j  .The  Log(error)  as  a  function  of  the 
number  of  iterations  for  the  same  problem 
done  three  different  times. 


Fig  13  The  Log(error)  as  a  function  of  the 
number  of  iterations  .The  same  problem 
done  twice  for  a  singular  matrix  A. 


have  infinite  solutions.  However,  the  BOC  can  be  used  to  solve  such  systems  [9].  1  he 

BOC  converges  much  faster  when  A  is  singular,  because  a  nonsingular  matrix  will  have  a 
unique  solution.  Due  to  the  infinite  solutions  that  a  singular  matrix  has,  the  BOC 
produces  different  solution  each  time  we  try  to  solve  the  same  problem  again.  Fig  13  shows 
the  BOC  convergence  for  a  singular  matrix. 

V.  Using  the  BOC  for  Phased  Arrays 

It  is  shown  the  previous  sections  that  the  BOC  is  capable  of  large  ill-conditioned 
linear  systems  of  equations  very  rapdily.  This  makes  the  BOC  a  unique  system  in 
processing  of  phased  array  data.  In  this  section  we  show  some  simulation  results  using  the 
BOC  for  the  weight  adaptation  of  a  phased  array. 

Two  simulation  experiments  are  presented  in  this  section.  In  the  first  experiment 
we  used  a  five  element  array,  and  assumed  the  directions  of  the  jammers  are  known.  In  the 
second  experiment  a  2  element  array  is  used  and  no  a  priori  information  is  assumed. 

In  Fig.  14  the  5  element  array  pattern  is  plotted  as  a  function  of  the  angle,  Fig. 
14(a)  shows  the  array  pattern  before  adaptation.  In  Fig.  14(b)  the  pattern  after  adaptation 
is  shown  for  a  jammer  at  45°.  The  array  pattern  after  adaptation  has  reformed  in  such  a 
way  to  null  the  jammer  signal.  In  Fig.  14(c)  four  jammers  are  considered  at  45°,  80  ,  120° 
and  150°,  the  array  pattern  is  again  reformed  to  null  all  the  jammers  signals  reception. 

In  Fig.  15  the  BOC  was  used  to  solve  the  adaptation  problem  assuming  no  a  prion 
information  about  t lie  interference  signals.  Fig.  15(a)  shows  the  two— element  array 
pattern  before  adaptation.  In  Fig.  15(b)  to  (d)  the  pattern  is  plotted  for  a  single  jammer 
placed  at  30°,  45°  and  60°.  respectively.  In  all  these  plots  the  array  adapted  to  cancel  the 
interference  signal  in  each  of  the  given  cases.  In  all  of  the  above  results  the  jammer  signals 
is  considered  to  tie  of  the  same  strength  ns  the  desired  signal,  and  the  convergence  of  the 
solution  obtained  in  less  than  five  iterations  Also  the  condition  number  of  the  R  is 
between  if/'  and  it. 


VI.  Conclusions 


Hybrid  optoelectronic  processors  were  demonstrated  to  be  used  for  processing  the 
adaptive  phased  arrays  data.  These  processors  were  shown  to  have  the  high  accuracy  of 
digital  processors  while  are  faster  in  solving  large  systems  of  linear  equations.  In  this 
particular  application  of  the  bimodal  optical  computer  considered  here,  it  has  been  shown 
that  the  HOC  is  suitable  for  the  adaptation  processes  of  the  weights.  The  interference 
canceling  was  achieved  even  with  the  presence  of  the  errors  in  the  analog  optical  processor. 
Experimental  results  reported  here  confirms  the  theoretical  predictions. 

For  future  work  the  I30C  can  be  implemented  to  do  real  time  adaptation  of  the 
weights.  Also  other  algorithm  adaptation  can  be  considered. 
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ABSTRACT 

A  novel  system  for  solving  systems  of  nonlinear  equations  is  proposed.  Two  diiferen: 
algorithms  are  introduced  A  speed  analysis  of  the  two  different  algorithms  is  presented  and 
compared  with  the  speed  of  their  digital  computer  counter  part.'.  A  great  advantage  m  speed  ].- 
shown  for  large  size  problems. 


1.  INTRODUCTION 

Systems  of  nonlinear  equations  arise  in  the  process  of  solving  many  physical  problem.'.  They 
are  a  very  important  class  of  mathematical  problems.  Iterative  methods  are  used  to  soive  such 
problems. 

In  this  paper  we  propose  a  new  method  for  solving  this  class  of  nonlinear  problems  using 
optical  processors.  In  Section  2  the  iterative  methods  used  in  solving  nonlinear  systems  m 
equation  is  reviewed.  In  Section  3  the  optical  implementation  is  proposed  using  two  different 
algorithms.  The  speed  analysis  of  the  two  algorithms  is  given  in  Section  4.  In  Section  5 
conclusions  and  final  remarks  are  drawn. 

2.  NEWTON's  METHOD 

Systems  of  linear  equations  are  given  as  follows 

-♦  -♦ 

Ax  =  b  A  ' 


-*  -t 

where  A  is  an  n  *  n  matrix,  x  and  b  are  n  x  1  vectors.  In  these  systems  A  and  b  ate  giv  u  ana  x 
the  solution  of  the  system  is  unknown. 

Nonlinear  systems  of  equations  can  be  represented  by 


or 


f(*)  =  o 
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where  fj's  are  nonlinear  functions  of  x 

One  of  the  methods  used  in  solving  for  x  in  the  nonlinear  system  of  equations  is  Newton's 
method.  For  a  single  nonlinear  equation,  an  initial  solution,  x®,  of  the  equation  is  assumed,  and 
the  (k-M)th  iteration  of  the  solution  is  given  by(1: 

xk+ !  =  *k  -  fk  <3> 

where 

fk=f(xk>'andfk  =  ^i)|  x  =  xk  (4) 

For  a  system  of  nonlinear  equations,  Eq.  (3)  can  be  rewritten  as 

*k+l  =  *k-(Jk)  1  ^k  ’  ^ 


where 

and  J  is  the  Jacobian  matrix. 

Let. 


(6) 


(7) 


men 


(8) 


Eo.  (SI  is  a  system  of  linear  equations  to  be  solved  for  2k,  which  is  the  correction  needed  for  the 
(k+I)th  solution  iteration.  The  algorithm  for  solving  the  system  of  nonlinear  equations  will  be  as 

loilows: 

i  I  Assume  a  solution  20. 

ii)  Compute  the  nxl  vector  ?k  and  the  nxn  matrix  Jk 

iii ,  Solve  the  linear  system  of  equations  Jk2k  =  ?k  for  2k 

i  v )  Compute  the  refined  solution  5tk+t  =  *k  - 

v )  Check  if  the  norm  ||?k*i~  ?kl|  <  c  stop,  otherwise  go  back  to  step  (ii).  c  is  the 

allowable  error. 

3.  OPTICAL  IMPLEMENTATION 

The  iterative  algorithm  introduced  in  Section  2  requires  0(n3)  number  of  operations  when 
used  with  conventional  digital  computer.  The  most  expensive  part  of  the  algorithm  is  step  (iii) 
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to  solve  a  system  of  linear 
equations.  In  previous 

publications^-4'  we  have 
proposed  and  analyzed  a  hybrid 
optoelectronic  processor,'  the 
Bimodal  Optical  Computer 
BOC,  capable  of  solving  linear 
systems  of  equations  accurately 
and  rapidly.  In  this  section  we 
modify  that  system  to  be  used 
to  solve  systems  of  nonlinear 
equations  as  shown  in  Fig.  1. 
We  propose  two  different 
algorithms,  the  first  utilizes  the 
use  of  the  analog  processor  to 
solve  the  system  of  equations 
(8)  approximately,  and  the 
second  to  use  the  BOC  to  solve 
the  system  of  equations  (8) 
exactly  (within  the  specified 
accuracy). 


Fig.l  Block  diagram  of  the  hybrid  optoelectronic  system. 


3.1  Hybrid  Analog  Optical  Processor 

In  this  system  we  use  the  optical  analog  processor  to  solve  Eq.  (S)  approximately.  For  this 
system  we  introduce  the  following  algorithm: 

-♦ 

a)  Use  the  digital  processor  to  guess  an  initial  solution  xG. 

b)  Use  the  digital  processor  to  compute  both  the  vector  ?k  and  the  matrix  Jk- 

c)  Use  the  optical  analog  processor  to  solve  the  system  Jk  2k  =  ?k  for  2£.  approximately, 

where  the  superscript  o's  denote  inaccuracies  in  optics  or  electronics. 

d  i  Use  the  digital  processor  to  read  2£  and  compute  the  refined  solution  3£k«i  =  -  *\- 

e)  Check  if  the  norm  ||?k«i  -  ?k||  <  « stop,  otherwise,  go  back  to  step  (b)  and  recycle. 

3.2  Hybrid  BOC  Processor 

In  this  system  the  BOC  is  used  to  solve  Eq. ( 8 )  exactly.  For  this  system  we  introduce  the 
following  algorithm: 

a)  Use  the  digital  processor  to  guess  an  initial  solution  x0. 

b)  Use  the  digital  processor  to  compute  both  and  the  matrix  J^. 

-*  *♦ 

c)  Use  the  BOC  to  solve  the  system  Jk  Ck  =  fk-  exactly  for  Ck 

d)  Use  the  digital  processor  to  read  2k  and  compute  the  refined  solution  5lko  —  2k  -  2k. 

e)  Check  if  the  norm  ||?k«i  -  ?k||  <  f  stop,  otherwise,  go  back  to  step  (hi  and  recycle. 
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4.  SPEED  ANALYSIS 


The  foliowing  speed  analysis  is  based  on  a  system  of  linear  equations  with  size.  n. 

4.1  Digital  Processor 

The  total  time  required,  T^,  to  solve  the  system  of  nonlinear  equations  using  a 
conventional  digital  processor  is  given  by 

3 

I  pq1  =  [~ 3  f2n(n+l)]T^jNjj  (-1 

wtiere 

T'di  =  the  time  needed  to  do  one  digital  operation  (e.g.,  a  multiplication). 

and 

ND  =  the  number  of  iterations  needed  for  the  solution  convergence. 

4.2  Hybrid  Analog  Optical  Processor 

The  total  time  required.  TQ. ,  to  solve  the  system  of  nonlinear  equations  using  the  processor 
introduced  in  Section  3.1  is  given  by 

Toa=  [n(n+2)Tu,  +  TA,|  NA  (10). 

where 

T A  j  thr  ..’  .  required  for  the  optical  analog  processor  to  solve  the  system  of  linear 

rquations  (8)  approximately, 
and 

N  ^  -  the  number  of  iterations  required  for  the  solution  convergence. 

4  3  Hybrid  HOC  Processor 

The  total  time  required,  Tq  ,  to  solve  tiie  system  of  nonlinear  equations  using  the  processor 
introduced  in  Section  3.2  is  given  by 

TOB  =  (2n(n+1)TDl  +  TAll  !BND  ' 1 1 1 

where 

1^  -  the  number  of  iteration  needed  for  the  BOC  to  solve  Eq.  (8)  to  the  specified  accuracy 
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4.1  Speed  Advantage 

It  is  of  great  interest  to  determine  what  is  the  break  even  point  for  the  optical  processor 
proposed  to  be  faster  than  the  digital  processors.  This  condition  is  defined  by 

tdt>toa  021 


tdt>tob 

From  Eqs.  (91  to  (11)  the  conditions  (12)  and  (13)  can  be  written  as 


n  (n/3+r 


fe) 


for  the  hybrid  analog  processor,  where 


And  for  the  hybrid  ROC  processor 


A  *  A  >  1 
n  t 


lA  ~  NA/% 


Where 


^Vn(n  +  l)(lB-i)  j  x  jhjl)  , , 


B  x  A.  >  1 
n  t 


n  (n/3+1) 


B  = 
n 


n3/3  -  2n(n+l)(IB-l) 


The  number  of  iterations.  1^  and  1^,  usually  are  in  the  range  of  1  to  ln  1  I  ■  r;  tins.  An  and  Bu 

are  problem  dependent,  and  are  much  larger  than  1  for  large  values  of  n.  On  me  other  nano.  At 
depends  on  the  speed  of  the  analog  processor  for  solving  a  system  of  linear  equation,  wtiicn  ean  no 
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in  the  ramie  of  /iscc.  But  since  the  matrix  Jk  need  to  be  updated  every  cycle,  writing  the  inatn.*. 
.Ik  on  the  SLM  becomes  the  bottleneck  of  the  processor  speed.  With  todays  technology  writing  ,1 
matrix  on  an  SLM  may  take  a  few  milliseconds.  So  A  is  much  iess  than  1.  In  Fig.  2(aj  and  i  Li 

the  Log  (An)  and  Log  (Bn)  are  plotted  in  terms  of  the  system  size,  n,  respectively.  I  he  ratio  AB  > 
1.  for  n  ~  10.  while  An  >  1  for  n  A  60  and  120.  for  IR  =  10  and  20  respectively.  For  the  At  ratio  in 

me  range  of  10'3,  we  can  have  a  speed  advantage  for  the  hybrid  analog  optical  processor  for  n  > 
50.  and  for  the  hybrid  BOC  processor  for  n  >  120. 


Fig. 2  Plot  of  log  of  the  ratio  (a)  Aq  of  Eq.(  19).  and 
(b)  Bn  of  Eq.(20),  in  terms  of  the  size  of  the  matrix  ,n. 


Again  this  ratio  At  depends  mainly  on  how  fast  we  can  write  a  matrix  on  the  SLM.  By  the 
introduction  of  faster  SLM's  the  speed  advantage  can  be  gained  for  smaller  values  of  n. 

5.  CONCLUSIONS 

;  wo  new  hybrid  opto  electronic  processors  are  introduced  for  solving  systems  of  nonlinear 
equation-  The  speed  of  the  two  processors  is  analyzed  and  compared  with  the  speed  of  digital 
processor/  It  is  shown  that  the  main  factor  of  the  speed  limitation  is  the  speed  the  SLM's  used  to 
write  the  matrix  on.  With  the  existing  SLM's  a  speed  advantage  can  be  gained  for  n>100. 
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ABSTRACT 

Hardware  and  software  design  of  the  Bimodal  Optical  Computer  (BOO  and  its 
implementations  are  presented.  Experimental  results  of  the  BOC  for  solving  a  system  of  linear 
equations  Ax  =  b  is  reported.  The  effect  of  calibration,  the  convergence  reliability  of  the  BOC. 
and  the  convergence  of  problems  with  singular  matrices  are  studied. 

1  .INTRODUCTION 

Analog  optical  systems  are  becoming  very  attractive  in  the  area  of  signal  processing  l .-■re use 
of  their  ability  to  process  in  parallel  two  dimensional  data  very  rapidly.  However,  analog  optica. 
systems  have  low  accuracy.  BOC  [1—4]  solves  this  low  accuracy  problem,  by  using  a  comoinatios. 
of  both  analog  optical  system  and  digital  processor. 

In  this  paper  we  present  experimental  results  using  BOC  for  solving  system.-  ul  linear 
equations.  In  Section  2  a  comparison  between  astigmatic  optics  and  waveguides  oared  algebra 
processors  is  presented.  The  hardware  and  the  software  design  of  BOC  is  in  Section  3.  Section  -5 
contains  the  experimental  results  of  the  BOC  for  solving  a  system  of  linear  equations,  i  he 
conclusions  are  in  Section  5. 

2.ASTIGMATIC  OPTICS  AND  WAVEGUIDES  BASED  ALGEBRA  PROCESSORS 

The  analog  optical  system  can  be  applied  in  many  applications.  This  paper  concern  rate?  on 
solving  a  system  of  linear  equations.  Goodman  [5]  has  introduced  an  astigmatic  processor  to 
perform  matrix  vector  multiplications,  which  can  also  be  used  in  a  system  of  linear  equations 
solver.  However,  the  main  problem  that  faces  the  arrangement  in  Fig.  1  is  aligning  the 
components,  to  insure  a  uniform  light  distribution  along  the  matrix  plane. 

Waveguides  can  be  used  to  build  optical  algebra  processors.  By  using  waveguide;,  the 
optical  system  can  be  made  compact,  and  its  alignment  will  be  much  easier  than  that  of  the 
astigmatic  system.  The  distribution  of  the  light  across  the  waveguide  is  plotted  in  Fig.  2.  which 
shows  that  the  light  is  almost  uniform  along  the  waveguide.  From  the  practical  standpoint 
waveguides  are  more  reliable  to  use  in  these  systems  than  the  astigmatic  optics. 
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3.THE  ROC  DESIGN  (HARDWARE  AND  SOFTWARE) 

3.1  BOC  HARDWARE  DESIGN 

The  BOC  hardware  has  three  main  parts  as  shown  in  Fig.  3.  The  optical  system,  the 
electronic  circuit,  and  the  digital  processor.  The  optical  system  consists  of  the  fully  parallel 
matrix-vector  multiplier.  Light  from  the  LED's  representing  the  x  components  are  spread 
vertically  by  planner  waveguides  onto  the  columns  of  the  matrix  mask.  The  transmitted  light  i? 
summed  row  wise  by  using  another  set  of  planner  waveguides  and  detected  by  photodiodes  which 
represent  the  output  vector  b. 

The  electronic  circuit  acts  as  a  feedback  loop  to  correct  for  the  input  light  of  the  LED's, 
until  a  solution  is  reached.  The  solution  x  will  then  be  read  and  stored  by  the  digital  processor. 
Fig.  4  shows  the  electronic  circuit  used  for  the  feedback  loop. 

The  A/D  and  D/A  conversion  from  and  to  the  electronic  circuit  are  performed  by  the  digital 
processor. 

3.2  BOC  SOFTWARE  DESIGN 

The  BOC  software  controls  the  Input/Output  operations.  Both  the  matrix  A  and  the  output 
vector  b  are  read  and  stored  by  the  digital  processor.  The  vector  b  is  then  converted  to  analog 
voltage  by  a  D/A  converter,  and  it  is  assigned  to  the  different  ports  of  the  electronic  circuit.  The 
analog  optical  processor  solves  for  an  approximate  solution  due  to  its  inaccuracy.  The  digital 
processor  reads  and  stores  the  approximate  solution,  x£  through  the  A/D  converter,  then  it 
calculates  the  residue  vector,  r,  as, 

r=b-Ax°=A(x-x°)=AAx  (1) 

Multiply  Eq.  (1)  by  a  scalar  s  to  make  use  of  the  whole  dynamic  range  of  the  system,  so 

F.q  1 ! )  becomes, 


sr=A(sAx)  (2) 

!f  the  residue  is  not  small  enough,  the  system  of  linear  Eq.(2)  will  be  solved  for  Ax  using  the 
analog  optical  processor  and, 

x*-x°+Ax  (3) 

A  new  residue  will  be  found  for  x*.  The  iteration  process  is  continued  by  solving  Eqs.(l) 
through  (3 1  until  a  satisfactory  solution  is  reached. 

4.EXPERIMENTAL  RESULTS 

In  this  section  we  present  the  experimental  results  for  solving  a  system  of  linear  equation? 
Ax=b  using  the  BOC, where  A,b,  and  x  are  all  positive. 

The  Log  of  the  error  and  that  of  the  residue  are  plotted  versus  the  number  of  iterations.  The 
error  and  the  residue  are  defined  as. 

Errors  j|(x-xk)||/||x||  (T 
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Residue=j|r^||  (5) 

Where  ||.||,  is  the  Enclidean  norm,  x  is  the  exact  solution,  xk  is  the  kth  iteration  solution, 
and  rk  is  the  kltl  iteration  residue. 

Since  we  are  dealing  onlv  with  positive  numbers  in  this  paper,  we  used  the  absolute  value  of 
r  to  solve  Eq.(2),  then  we  set: 

x(n+l)=x(n)+Ax  (6) 

when  all  the  components  of  r  are  positive.  And 

xl“+1WnFlx  (') 

if  all  the  components  of  r  are  negative.  We  reject  the  iteration  when  the  components  of  r  have 
different  signs  and  take  the  previous  one.  By  rejecting  some  iterations  we  are  actually  rejecting 
some  corrections.  This  procedure  slows  down  the  convergence  process. 

In  all  the  experiments  performed,  the  iteration  process  is  stopped  when  a  16  bit  accuracy  is 
reached.  Fig.  5  shows  that  the  BOC  started  with  almost  30%  error  and  it  needed  6  iterations  to 
converge  to  16  bit  accuracy.  In  Fig.  6(a)  BOC  started  with  almost  110%  error,  and  the  number  of 
iterations  needed  was  21.  Fig.  6(b)  shows  the  Log  of  the  residue  as  a  function  of  the  number  of 
iterations.  The  fluctuations  depicted  by  Figs.  6(a)  and  (b)  is  due  to  the  rejection  method  used  in 
the  experiments. 

4.1  EFFECT  OF  CALIBRATION 

The  analog  optical  system  error  is  a  major  factor  in  the  rate  of  convergence  of  the  BOC.  If 
that  error  is  reduced,  then  the  convergence  is  much  faster.  In  order  to  illustrate  this,  the  Seme 
problem  has  been  solved  twice  with  two  different  accuracies  of  the  optical  system.  The  analog 
optical  system's  error  in  the  first  time  was  50%,  and  it  was  30%  in  the  second  time.  Twer,  y  one 
iterations  were  needed  by  BOC  to  converge  to  the  16  bit  accuracy  for  the  first  case  for  trie 
second  case  the  number  of  iterations  was  reduced  to  12.  These  results  are  plotted  in  Fig 

4.2  RELIABILITY  OF  THE  SYSTEM 

System  reliability  for  convergence  have  been  tested  and  verified  by  solving  the  same  problem 
several  times,  under  different  conditions.  Results  show  that  when  the  BOC  is  used.  *o  solve  a 
problem  several  times,  the  convergence  rate  will  not  be  exactly  the  same  for  all  the  cases. 
However,  the  number  of  iterations  needed  by  the  BOC  to  converge  to  a  certain  accuracv  is  almost 
the  same.  Fig.  8  shows  three  different  paths  of  convergence  for  the  same  problem,  i  no  BOC 
needed  13  iterations  in  the  first  run,  14  iterations  in  the  second,  and  11  in  the  third. 

4.3  SOLUTION  CONVERGENCE  FOR  THE  SINGULAR  MATRIX  SYSTEM 

Solving  a  system  of  linear  equations  with  a  singular  matrix  A  is  one  cf  the  problem?  that 
cannot  be  solved  using  conventional  digital  computer  techniques.  Singular  matrices  have  a 
condition  number  equal  to  infinity,  so  their  inverse  does  not  exist,  also  they  have  infinite  number 
of  solutions.  However,  the  BOC  can  be  used  to  solve  such  systems  [(>'.  The  BOC  converges  much 
faster  when  A  is  singular.  because  a  nonsingular  matrix  will  have  a 
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unique  solution.  Due  to  the  infinite  solutions  that  a  singular  matrix  has.  the  BOC  produces 
different  solution  each  time  we  trv  to  solve  the  same  problem  again.  Fig.  9  shows  the  BOC 
convergence  for  a  singular  matrix,  r 


5.CONCLUSIONS 

The  BOC  system  was  built  and  experimentally  tested.  The  experimental  results  show  great 
reliability  of  the  processor  in  solving  systems  of  linear  equations.  Overall  16  bit  accuracy  of  the 
hybrid  system  was  achieved  with  an  analog  optical  system  of  30%  to  50%  error.  Higher  accuracies 
of  the  solution  can  be  obtained  by  increasing  the  number  of  iterations.  The  BOC  also 
demonstrated  to  solve  systems  of  linear  equations  with  singular  matrices. 

We  are  considering  in  future  work,  bipolar  numbers,  complex  numbers,  and  using  SLM  for 
the  matrix  mask. 
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Fig  1.  An  astigmatic  processor  for  solving 
system  of  linear  equations. 
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Fig  2.  Light  distribution  acrof;,  the 
waveguide. 


•  nxcTRONic  cincvrr 


Fig  3.  1  he  Bimodal  Optical  Computer. 


Fig  4.  circuit  diagram  of  the  feedback 
loop  for  the  optical  system. 
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Fig  5.  The  Log(error)  as  a  function  of 
the  number  of  iterations.  The  BOC  started 
with  30%  error. 


Fig  8.  The  Log(error)  as  a  function  of  the 
number  of  iterations  for  the  same  problem 
done  three  different  times. 
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Fig  6(a).  The  Log(error)  as  a  function  of 
the  number  of  iterations.  The  BOC  started 
with  100%  error. 


Fig  7.  The  Log(error)  as  a  function  of  the 
number  of  iterations  for  the  same  problem, 
but  with  two  different  accuracies  of  the 
optical  system. 


Fig  9.  The  Log(error)  as  a  function  of  the 
number  of  iterations  .The  same  problem 
done  twice  for  a  singular  matrix  A. 
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ABSTRACT 

A  set  of  ill-posed  aigebra  problems  are  considered  for  solving  by  > !;<•  B;  modal  (Erao. 
Computer.  BOC.  The  BOC  algorithm  was  shown  to  be  capable  of  solving  ihis  class  o:  *i 

problems.  Three  different  methods  of  generating  the  error  matrix  are  compared  m  term-'  o: 
convergence  of  the  solution.  Some  applications  for  the  methods  are  introduced. 

1.  INTRODUCTION 

Optical  linear  algebra  processors1  ■  are  introduced  to  provide  a  fast  and  reliable  !..• 
solving  linear  algebra  problems.  We  have  shown  in  previous  publications  -'  that  tin*  BO’” 
capable  of  solving  systems  of  linear  equations  with  singular  matrices.  The  BOC  alsoritnr.  v  '  !. 
we  have  introduced  can  also  be  used  on  conventional  digital  computers.  In  this  paper  we  ptesot-.t 
three  different  methods  in  generating  the  error  matrix  used  in  the  algorithm.  A  commit :<c 
between  the  different  methods  is  shown  in  terms  of  the  solution  convergence. 

In  Section  2  we  review  the  BOC  algorithm  for  solving  a  system  of  linear  equation." 
Different  schemes  in  generating  the  error  matrix  is  given  in  Section  3.  Application  u:  moth 

to  semi  linear  problems  is  shown  in  Section  1.  In  Section  5  conclusions  and  final  remarks  2:ve:. 

2.  THE  BOC  ALGORITHM 

In  this  section  we  review  the  BOC  algorithm  to  solve  a  system  of  linear  equation.-  ! 
algorithm  originally  was  developed  for  the  Bi modal  Optical  Computer. BOO  wnid.  :  a  ;:\m  ■ 
system,  combining  both  analog  optics  and  digital  electronics  to  achieve  accurate  -lint tom  vc:\ 
rapidly  for  systems  of  linear  equations.  The  algorithm  can  also  be  used  on  a  diciia.  con. tvitei 
solve  systems  with  very  iil-conditioned  or  singular  matrices 

the  problem  we  are  interested  to  solve  is  the  basic  system  of  linear  equatio.  • 


Ax  =  b  . 


where  A  is  an  n  *  n  matrix,  x  ami  b  are  n  *  1  vectors.  A  aim  b  are  are  Known.  aim  x  w 
solution  of  the  system  we  need  to  determine.  The  matrix  A  can  be  eitimi  an  i i i  — « < u:- : : t : ; o •<  1  o: 
-jugular. 

The  BOC  algorithm  is  as  follows: 
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U'i  Generate  the  matrix  A0,  which  is  given  by 

A0  =  A  +  E  . 


whore  E  is  the  error  matrix. 


Pi  • 

•  tin 

E  = 

til  ' 

*  l  IQ 

t  n  1  • 

‘  Cqq 

L 

and  is  generated  using  a  Gaussian  random  number  generator,  with  zero  mean  and  standard 

deviation,  ac. 

E 

(b''.  Solve  the  system 

A  x  -  b  ,  14) 


icj  Compute  the  residue 


r  =  b-  Ax-AAx. 


(5) 


!d)  Solve  the  system 


-  A°  sAi° 


Ax'  .  where  s  is  a  constant  used  to  utilize  the  dynamic  range  of  the  system,  such  that 


(6) 


e  Compute  the  refined  solution 


P  =  i°  +  At0  (S' 

* 

C  if  ||  r  jj  <  e.  the  allowable  error,  stop,  otherwise,  go  to  step  (c)  and  reiterate. 

In  this  algorithm  the  matrix  A°  is  a  well  conditioned  matrix* 4 ' .  If  A  is  a  singular  matrix 
then  it.  has  at  least  one  of  its  eigenvalues  is  equal  to  zero.  By  adding  the  error  matrix  E  to  A  then 
‘he  zero  eigenvalues  of  A  will  be  shifted  ami  the  matrix  A0  will  have  no  zero  eigenvalues.  It  is 
also  worth  mentioning  here  that  step  (a)  is  not  necessary  when  using  an  analog  optical  processor, 
because  it  is  done  naturally. 

3.  GENERATING  THE  ERROR  MATRIX 

In  this  section  w.-  introduce  other  ways  of  generating  the  error  matrix.  E.  given  in  Eq.  :  J  ■ . 
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Also  we  wiii  compare  these  different  methods  in  terms  of  the  solution  convergence  of  svstem  of 
linear  equations  with  a  singular  matrix. 

The  error  matrix.  E.  can  be  generated  as  given  by  Eq.  (3).  Also  some  of  the  other  means  it 
can  be  generated  are 


or 


r<  o 


or 


10 


fi 

0 

ro 

1 

•  1 

+  c0  = 

•  0 

[0 

■  1 

£ 

1 

■  0 

In  Eq.  (9)  the  diagonal  matrix  is  generated  using  Gaussian  statistics. 

In  Fig.  1  the  number  of  iterations  required  for  a  16  bit  solution  of  a  5  *  5  singular  svstem 
plotted  in  terms  of  the  logarithm  of  the  standard  deviation.  <7p,  of  matrix  E  and  Ej  or  the 

logarithm  of  c  of  matrix  E^.  In  Fig.  1(a)  a  5  *  5  matrix  with  rank  =  4  is  used,  while  in  Fur.  hi, 
the  rank  of  the  matrix  is  only  1. 


Fig.  1  Number  of  iterations  required  for  16  bit  resolution  is  plotted 
in  terms  of  the  logarithm  of  or  c  for  a  5x5  matrix  with 

fai  rank=4.  and  t'b)  rank  =  l. 

In  Fig.  F  a)  the  number  of  iterations  using  E  and  Fl|  is  5  for  an  ernumb'T  of  iteration  starts 
;o  increase.  While  the  number  of  iterations  using  Eg  is  only  4  for  an  error,  f.  up  to  lO'-c.  For  rank 
I  matrix.  Fig.  Kb),  the  number  of  iterations  is  3  or  4  for  an  error  up  to  id'  '.  When  the  erto; 
roar  lies  100'/  the  number  of  iterations  is  between  15  and  25.  with  t  h  *  *  method  using  Eo  has  the 
best  performance. 
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In  Fin  2  tin1  number  of  iterations  required  for  a  16  bit  resolution  solution  is  plotted  m  term 
;  Log(Ci/C-„>).  for  two  values  of  C->.  I' he  number  of  iterations  is  5  up  to  an  error  of  1000a 


except  a  divcigence  at  Ct  =  C?.  whicii 
mai.o  the  error  matrix  itself  singular. 

From  ti.e  resuits  shown  in  Figs.  1 
and  the  different  means  of  generating 
aw  error  matrix  are  quite  similar  in  terms 
wf  the  convergence  rate.  It  is  clear  that 
.my  of  the  above  methods  is  capable  of 
■oivmg  tins  ciass  of  problems.  I  sins  an 
1 1  it  ic a!  processor  will  limit  us  to  use  E  for 
•no  error  matrix.  When  the  matrix  .A  .is 
w:  it  ten  in  the  analog  processor  because  of 
die  inherent  inaccuracies  of  the  processor 
A°  will  be  recorded  instead  of  A  But 
when  using  tins  algorithm  on  a  digital 
computer,  we  are  not  limited  to  any  one 
of  the  a  Dove  methods.  In  this  case  using 
E>  for  generating  the  error  matrix  will 
inquire  the  least  number  of  operations. 


F o  2  Number  of  iterations  required  for  16  bit 
resolution  in  terms  of  Logf C1/C2)  for  a 
5x5  matrix  with  rank=4. 


4.  SEMI  LINEAR  PROBLEMS 

This  ciass  of  problems  lias  a  very  wide  range  of  applications:  such  as  optimization  problems. 
:;;iy,r  programming,  image  restoration,  etc.  In  semi  linear  problems  the  system  of  it  near  equation 

iv.-u  in  Eq.  i ! !  need  to  be  solved  for  *  with  the  following  constraints 

Xj  >  0  ,  i  =  1 , -  -  - , n  (12) 


bj  >  0  .  (13) 

Usually,  the  system  of  linear  equations  considered  in  the  semi  linear  problems  is  overdetermined, 
underdetermined  or  with  matrix  A  singular.  In  all  these  cases  the  system  does  not  have  a  unique 
solution.  So  in  the  iterative  BOC  algorithm  any  solution  that  does  not  satisfy  Eq.  (12)  will  be 
icpiaceu  zero,  until  the  system  converges  to  a  satisfactory  solution. 


5.  CONCLUSIONS 

Fiie  BOC  algorithm  is  introduced  and  applied  to  solve  a  set  of  ill-posed  algebra  problems. 
Tiw  convergence  of  the  solution  was  studied  for  the  different  methods  used  to  generate  the  error 
niatiix  The  algorithm  can  be  used  both  optical  and  digital  computers. 
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ABSTRACT 


The  Bimodal  Optical  Computer  (BOC)  is  considered  for  Adaptive  Phased  Array  Radar  (APAP.)  data  processing.  The 
effect  of  the  errors  in  the  BOC  on  the  optimum  weight  calculations  for  the  interference  canceling  are  studied.  Computer 
simulations  for  five  and  nine  element  APARs  are  presented. 

I.  INTRODUCTION 

The  adaptive  phased  array  radar  systems  provide  the  means  for  suppressing  unwanted 
interference  signals.  This  is  achieved  bv  nulling  the  antenna  pattern  at  tne  direction  of  the 
jammers.[l-5]  Many  algorithms  are  introtfuced  for  the  adaptation  process  and  these  are  reviewed 
by  Monzingo  and  Miller.[2]  These  algorithms  have  large  computational  complexity.  This  makes 
the  adaptation  process  slow,  especially  for  large  size  arrays. 

We  have  introduced  the  oimodal  optical  computer  (BOC)  for  solving  linear  algebra 
problems.(6-7]  BOC  is  capable  of  solving  tnis  class  of  problems  with  both  high  speed  and  accuracy. 
As  an  optical  hybrid  system,  the  BOC  combines  between  the  speed  and  parallelism  of  optical 
analog  processors  and  the  accuracy  of  the  digital  electronics.  In  previous  paper  we  have  suggested 
the  implementation  of  the  BOC  for  processing  of  the  APAR’s  data.[8] 
in  this  paper  we  study  the  effect  of  the  errors  involved  with  the  BOC  components  on  the 
interference  canceling  ability  of  the  APAR.  Also  a  comparison  between  the  antenna  pattern  for  an 
AFAR  computed  using  the  BOC  with  and  without  optical  errors  is  presented. 

1 . 1  Adaptive  Phased  Array  Radar 

The  adaptive  phased  array  radar  system  (APAR)  block  diagram  is  shown  in  Fig.l.  The  received 
signal  vector  x  is  given  by 

X  =  X<j  +  Xj  +Xn,  (1) 

where  x^  =  the  desired  signal, 

x.  =  the  interference  signal,  and 
xn  =  the  thermal  noise. 

The  received  signal  Xj’s  are  then  multiplied  by  the  corresponding  weights  wj’s.  and  summed  to  give 
the  output  signal  s(t).  The  output  signal  is  then  compared  with  the  reference  signal  r(t)  to  give  the 
error  signal  <  (t). 

<  (t)  =  r(t)  -  wT  x,  (2) 

where  denotes  the  transpose  of  the  vector. 

The  adaptation  problem  for  the  weights  can  be  formulated  in  such  a  way  that  the  least  mean 
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square  error 


E{e2(t)}  =  E{[r(t)-wTx]2}  (3) 

is  minimized,  where  E{  .  }  is  the  ensemble  average.  In  the  static  case  the  least  mean  square  error 
(LMS)  can  be  set  to  zero.  This  can  be  achieved  by  solving  the  system  of  linear  equations  given  by 
[2], 

^xx  w  =  rxd  •  W 

Where  R^  =  E{x  x^  },  and  rxcj  =  Ejx^ r  }.  So  the  problem  for  computing  the  optimum  weights 
w  is  reduced  to  the  solution  of  a  system  of  linear  equations  given  by  Eq.  (4).  The  covariance  matrix 
is  symmetrical  and  usually  either  singular  or  ill-conditioned.  Solving  such  system  of  equations 
with  limited  accuracy  processor  will  require  a  very  large  number  of  iterations,  which  in  turn 
increases  the  computational  time. 


In  this  subsection  we  review  briefly  the  BOC  algorithm  which  is  used  to  solve  the  adaptation 
iroblem  given  by  Eq.(4).  Let  A  =  R^,  x  =  w,  and  b  =  rxCj.  Substituting  in  Eq.(4)  will  reduce  it  to  the 
ollowing 


Ax  =  b  .  (5) 

Which  the  standard  format  for  a  system  of  linear  equations.  The  BOC  block  diagram  is  shown  in 
Fig.2.  It  is  an  optical-hybrid  system.  The  optical  analog  processor  is  used  in  solving  the  system  of 
linear  equations  of  Eq.(5)  and  a  digital  processor  to  refine  the  accuracy  of  the  solution.  The 
algorithm  used  in  solving  the  system  of  equations  is  given  as  follows: 

(a)  Solve  Ax  =  b  using  the  analog  optical  processor  to  get  x0. 

(b)  With  a  dedicated  digital  electronics  processor,  read  x0  and  evaluate  the  residue 

ro  =  b-Ax0  =  Aax0  (6) 


(c)  Normalize  tv,  to  use  the  dynamic  range 

(d)  Solve  opticahy  the  system  of  equations 


where 


A z  =  srQ  , 


z  =  s  (ax0), 

s  =  radix  used  in  normalizing  r. 
(e)  Evaluate  electronically 


(7) 

(8) 


xi  =  *0  +  a*0  ,  (9) 

and 

r,  =  b  -  Axj .  (10) 

(f)  If  ||  rt  ||  is  small  enough  stop.  Otherwise  go  back  to  step  (c)  and  recycle. 

There  are  a  number  of  error  sources  in  the  BOC.  There  are  reading  and  writing  errors  in  x  and  b 
and  errors  in  writing  the  matrix  A.  Since  the  processor  is  analog,  thus  inherently  it  has  a  low 
accuracy  which  is  determined  by  the  dynamic  range  of  all  of  these  optoelectronic  devices.  These 
errors  will  affect  the  rate  of  convergence  of  the  solution,  which  also  depends  on  the  condition 
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number  of  the  matrix  A.  In  the  following  section  we  study  the  effect  of  these  errors  on  the 
performance  of  the  APAR. 

II.  COMPUTER  SIMULATION  RESULTS 

We  have  conducted  a  set  of  numerical  experiments  to  simulate  the  BOC  in  calculating  the 
optimum  weights  for  the  APAR.  The  simulation  experiments  are  considered  with  and  without 
errors  encountered  in  the  BOC.  This  is  done  to  study  the  effect  of  the  optical  errors  in  calculating 
the  weights  and  the  antenna  pattern.  In  the  following  results  we  have  considered  the  errors  to  be 
Gaussian  distributed  with  a  standard  deviation  of  1%  in  the  matrix  A  and  the  vectors  x  and  b. 


2.1  Five  Element  Array 

In  the  first  simulation  experiment  we  considered  a  five  element  adaptive  phased  array  antenna. 
The  antenna  elements  are  spaced  by  half  a  wavelength.  In  all  the  cases  reported  the  strength  of  the 
desired  signal  strength  is  considered  to  be  equal  to  that  of  the  jammer.  In  Fig.3  the  antenna  pattern 
of  the  APAR  is  plotted  in  terms  of  the  azimuth  angle.  The  desired  signal  is  at  zero  degrees.  In 
Fig.3(a)  the  antenna  pattern  is  plotted  for  a  five  element  array  with  one  jammer  at  30° .  The 
pattern  was  plotted  first  assuming  no  errors  in  the  BOC  (the  continuous  curve).  The  other  three 
curves  are  for  the  BOC  with  the  error  specified  above  for  1  iteration  (dashed  line),  2  iterations  (*) 
and  5  iterations  (x).  The  antenna  pattern  nulls  at  30° .  which  is  the  jammer  direction.  The  pattern 
calculated  witn  only  one  iteration  nulls  also  at  the  30°  but  not  as  deep  as  the  exact  one,  but  the  2 
and  5  iterations  patterns  null  as  deep  as  the  no  error  curve.  In  Fig.3(b)  the  same  APAR  is 
considered  with  a  thermal  noise  of  standard  deviation  equal  to  0.1  of  the  desired  signal  strength. 
The  continuous  curve  represents  the  pattern  obtained  by  the  BOC  with  optical  errors,  while  the 
dashed  line  is  the  BOC  with  the  errors  mentioned  above.  The  two  patterns  are  very  close  except 
that  the  dashed  one  is  not  as  deep  as  the  continuous  one,  with  one  or  two  extra  iterations  we  can 
expect  the  same  null.  Figs.3(c)  and  (d)  the  APaR  is  considered  with  2  jammers  at  30°  and  50° . 
Again  the  patterns  in  Fig.3(c)  is  without  thermal  noise,  the  nulls  are  very  deep  and  similar  with  no 
errors  or  with  a  practical  BOC  with  few  iterations.  In  Fig3(u;  thermal  noise  is  again  considered. 


2.2  Nine  Element  Array 

In  this  set  of  simulation  experiments  we  consider  9  element  APAR  with  a  number  of  jammers 
with  equal  strengths  to  that  of  the  desired  signal  which  is  at  0° .  In  Fig.4(a)  a  jammer  at  50°  is 
considered  and  no  thermal  noise.  The  pattern  is  plotted  using  the  BOC  with  no  errors  and  the 
BOC  with  errors  for  one,  and  three  iterations.  From  the  pattern  it  is  clear  that  the  pattern  plotted 
using  three  iterations  has  a  null  as  deep  as  that  of  the  system  with  no  errors.  In  Fig.4(b)  a  thermal 
noise  of  standard  deviation  equal  to  0.1  of  the  desired  signal  strength  is  considered  with  a  jammer 
again  at  50° .  The  antenna  pattern  obtained  using  the  BOC  with  no  errors  and  with  errors  and  one 
iteration  are  both  nulling  the  jammer  with  almost  the  same  depth.  In  Fig.4(c)  five  jammers  are 
considered  at  10° ,  20° ,  30° ,  40° ,  and  50° ,  and  thermal  with  strength  equal  to  that  of  the  desired 
signal.  Again  the  pattern  and  the  depth  of  the  nulls  are  very  similar  for  BOC  with  no  errors  or  with 
errors  after  few  iterations. 


III.  CONCLUSIONS 

The  bimodai  optical  computer  because  of  its  hybrid  nature  is  superior  in  speed  to  that  of  the 
digital  computer  in  solving  a  system  of  linear  equations,  especially  for  large  size  systems. [6]  It  is 
shown  here  that  the  BOC  is  capable  of  determining  the  optimum  weights  for  an  APAR  in  few 
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iterations  even  when  their  exists  large  inaccuracies  in  the  optoelectronic  devices  used.  Since  only  a 

very  small  of  number  iterations  are  needed  then  the  weights  can  be  updated  very  fast  which  is  a 

must  for  such  system. 
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Fig.  1  The  adaptive  phased  array  radar 
system. 


Fig. 2  The  bimodal  optical  computer  block 
diagram. 
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Fig.3  The  antenna  pattern  for  a  five  element  array  is  plotted  in  terms  of  the  azimuth  angle  for 
(a)  one  jammer  at  30°  and  no  thermal  noise,  (b)  one  jammer  at  30°  and  thermal  noise, 
( c)  two  jammers  at  30°  and  50°  and  no  thermal  noise,  and  (d)  as  in  (c)  with  thermal 
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Fig.4  The  antenna  pattern  for  a  nine  element  array  is  plotted  for  (a)  one  jammer  at  50°  with 
no  thermal  noise,  (b)  same  as  in  (c)  with  thermal  noise,  and  (c)  five  jammers  at  10°, 
20° ,  30“ ,  40“ ,  and  50“  with  thermal  noise. 
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ABSTRACT 

Tne  use  ot  the  himodal  optical  computer  1BOC1  in  determining  the 
ueignts  lor  an  aaapuue  pnasea  array  raaar  is  introduced  Intenerence 
canceling  is  presented  lor  mo  cases  ll  i  assuming  me  direction  ot  me 
tummer  is  known  and  ill  assumtnt  no  a  prion  miormalion  The  effect  ot 
me  rammers  on  me  array  pattern  is  snown  lor  up  to  lour  tammers. 

1.  INTRODUCTION 

The  senstuvttv  of  a  signal-receiving  antenna  array  svstem  to 
interfering  noise  sources  can  be  reduced  bv  suitable  processing 
of  the  outputs  of  the  individual  arTav  elements.  The  processing 
of  the  output  of  the  arrav  svstem  acts  as  an  adaptive  nltenng 
svstem  [1-4]  The  adapuve  phase  arrav  radar  systems  provide 
uie  means  of  suppressing  unwanted  interference  signals.  This 
is  achieved  ov  nulling  the  array  pattern  in  the  direction  of 
the  tammers  Manv  algorithms  have  been  introduced  for  the 
adapiauon  process  and  thev  are  reviewed  bv  Monungo  and 
■Miller  (2) 

In  this  paper  we  present  a  new  technique  for  determining 
the  weights  for  the  adapuve  arrav  using  the  bimodal  opucal 
computer  (BOC)  [5-7]  The  bimodal  opucal  computer  is  capa¬ 
ble  of  solving  svstems  of  linear  equation  vcrv  rapidly  with 
high  accuracv  In  the  adapiauon  process  we  reduce  the  prob¬ 
lem  to  a  svstem  of  linear  equauons.  which  in  tuna  is  solved 
using  the  BOC 

In  Secuon  2  we  review  the  basic  theory  of  adapuve  phased 
arrav  radars  The  bimodal  opucal  computer  algonthm  for 
solving  uie  adapiauon  problem  is  presented  in  Secuon  3 
Computer  simulation  results  are  given  m  Section  4  Conclu¬ 
sions  and  final  remarks  are  given  in  Section  5 

2  ADAPTIVE  PHASED  ARRAYS 

in  adaptive  phased  arrav  radars  the  incoming  signal  is  de¬ 
tected  bv  an  arrav  of  sensors  The  detected  signal  is  a  combi¬ 
nation  of  the  target  signal  plus  interference  and  noise  signals 
The  svstem  is  adiusted  in  such  a  wav  to  supprrss  the  inter¬ 
ference  signal  reception  without  atfecung  the  desired  signal 
In  this  section  we  consider  the  two  general  cases  of  inter¬ 
ference  canceung:  (li  bv  assuming  that  the  interference  signal 
direcuon  is  known  and  (2|  bv  assuming  no  a  prion  informa¬ 
tion  is  gnown  about  the  interference  signal 

Imerierence  Signal  Direction  is  Known  When  the  inter¬ 
ference  signal  direcuon  is  known  the  weights  >*■  of  the  arras 
can  be  cnosen  to  suppress  the  interference  signal  Let  the 
svstem  snosvn  in  f-igure  1(a)  be  used  to  demonstrate  this 
adapiauon  technique  The  output  signal  of  the  arrrav  it  n  is 
given  bv  i  1  \ 

1 1  ;  l  -  P\  l  u  -  w,  i  sin  uy.i  -  i»  »  »,  i  sint  us  t  -  "  -  ■  -  I  ] 

-  i  i  vini  ui.t  -  H  i  -  *.  sinl  w.  i  *  -  ■  -  i 

-  x.  sint  u  i  ~  (  i  -  x,  sint  ui.i  -  t*  -  ‘ "  l  j  .  (If 


wnere 

P  "  the  pilot  signal. 

/  *  the  imerierence  signal 
ft  -  me  pnase  snifi.  (  ; , 

Ira 

b  -  - sin  v 

\ 

To  cancel  me  imerierence  signal  ana  10  make  me  signal 
st  n  equal  to  the  pilot  signal,  we  need  to  solve  the  following 
svsiem  of  unear  equauons  for  me  weighis  x 

x-  x  »,  ”  0. 

(  3) 

ix  -  x  icos  0  —  I  x.  -  > sin  b  -  u. 

(  w.  m  x,  icos  0  -  (  X  -  x.  I  Sin#  “  0 

The  size  of  this  svsiem  of  linear  equauons  depends  on  the 
number  of  sensors  tn  tne  arrav  The  number  of  tammers  can 
make  the  svstem  under  or  overdetermined.  both  of  wrucn  are 
time  consuming  algebra  problems 

3.2.  .Vij  A  Prion  Information  is  Known  This  is  ihe  most 
general  case  wnere  we  assume  no  information  about  tammers 
The  svstem  used  in  this  case  is  shown  in  Figure  lib)  Each  of 
the  n  sensors  receives  a  signal  nil  that  is  in  turn  rr.uiupbed 
bv  a  variable  weighi  x  The  output  signal  s <  r  i  is  compared 
with  the  desired  signal  du)  and  their  difference,  the  error 
signal  till,  is  used  to  determine  the  value  of  x  The  output  of 
the  arrav  is 

J<r|  “  L  *,(/>*  (4| 

1 

or 

sfr)  -  w\.  (5) 

where 

'»t  ")  i  .i.frli 

w  -  »  |  and  \  *  i  till  i  (61 

For  digital  sampled  data 

if  / 1  -  w  x<  /  I  (  ') 

and 

<1  i )  -  d{  j  I  -  r(  j )  -  i/i  / 1  -  w  \(  i  |  ( ft) 

The  optimum  value  of  the  weights  «  is  the  one  that  reduces 
ti  j  i  to  zero  or  at  least  minimizes  it 

For  A  samples  of  data  me  optimum  weignts  satisfy  the 
following  set  of  svstems  of  linear  eauauons 

»  xl  1 1  -  hi  1 1 

( V) 

"  till  -  dill 
x  x(  V  |  -  dl  V  i 
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SIGNAL 


OUTPUT  y  (t ) 


Figure  1  Basic  aaaouvc  arrav  svsiem  witn  lai  sianai  ana  noise  directions  xnown  and  ibi  no  a  enon  imormanon  assumed 


The  ft  sci'.  of  eauauons  nave  n  unknowns,  and  usuaih 
\  »  n.  and  arc  inconsistent  and  ovcrsDecincd  The  ootirruza 
uon  problem  can  oc  rewritten  as 

w„r.  -  rtJ.  (10) 

where 

R  -  E |xx;  ;  (11) 

and 

r  .  -  Elxd  ;  (12) 

The  matnx  R  is  called  the  covariance  matrix  unere  El  i  is 
the  ensemble  averaec 

Manv  aJaonthms  are  introduced  I2|  lo  solve  lor  the  weiehts 
in  Eq  (10)  Some  ol  tne  pooular  aJaonthms  are  ine  least  mean 
square  (LMSi  and  tne  direct  matnx  inversion  (DMI) 

0»e  wall  bneriv  mention  the  DMi  aisoninm  since  n  leads  to 
the  aJtonthm  introduced  in  inis  paper  Equation  dOi  cannoi 
ne  determined  exacilv  usina  a  unuted  numocr  oi  samoies  o! 


the  input  data  For  practical  consideration  a  small  numocr  o: 
samples  is  detected  to  ne  used  in  dctermminE  w  The  esti¬ 
mated  vaiue  of  Eq.  ilui  can  oe  Eivcn  os 

“  =  R,,:r,_  .  ( 1 3  I 

where  ft.,  is  the  sample  covanance  matnx  and  r  ,  is  the 
sample  cross-correlation  vector  mat  are  erven  bs 

1  k 

R  «  -  _  \l  /  lx  ’(  /  j  114) 

■*v  _ 

and 


K  is  me  numoer  oi  simoies  Die  DMI  alEoninm  determines 
the  inverse  ot  the  sainme  covanance  matnx  K  and  then 
from  Eq  1 1  ’ )  evaluates  w 
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3.  THE  8IMODAL  OPTICAL  COMPUTER  ALGORITHM 

Convenience  of  either  the  LMS  or  the  DMI  aigontnms  de¬ 
pends  on  a  numner  of  factors,  tne  most  important  Peing  the 
conaition  numoer  of  the  malax  R,.  If  the  matrix  R.,  is 
til-conditioned  or  singular,  it  either  converges  verv  slowlv  or 
the  inverse  does  not  exist.  respectively,  in  such  cases  other 
methods  might  be  used,  but  thev  are  lengthiv  and  time  con¬ 
suming.  so  thev  are  not  suitable  for  a  svstem  where  ume  is  a 
verv  crucial  element 

We  nave  snown  in  previous  publications  16.  7)  that  tne 
oimoda)  optica)  computer  is  capable  of  solving  such  problems, 
wnere  the  svstem  of  equations  is  ill-condiuoned.  singular. 
.'versDccihed.  or  undersoecined  The  BOC  is  a  hvbnd  svstem 
■’v  nature,  see  Figure  2.  It  uses  analog  optics  to  solve  the 
problem  approximately  but  rapidlv  and  it  utilizes  the  digital 
electronics  to  renne  the  solution,  in  an  uerauve  seneme 

The  adaptation  problem  for  tne  weights  w  introduced  in 
Section  2.  con  oe  rewntien  in  tne  following  torm.  from  Eu 

R  w  «  r .  ( 1  b  l 

wfnen  can  oe  written  as 

■lx  -  b.  ( 1 7 1 

w  here 

4  -  R  . 

\  "  » 
b  =■  ( 

Equation  tin.  a  svstem  oi  iincar  cauations  mat  can  be 
volveo  using  me  nimodal  omicai  computer  Among  tne  ao 
vantages  ot  using  me  BOC  over  me  conventional  teenmauc's 
■ire  speed,  .evpeciallv  for  urge  size  arravvi.  convergence  of  me 
solution  tor  difficult  prooicms.  and  lii  conditioned  singular 
wviemv  union  is  me  cave  lor  movi  of  me  aoarioe  arras  raojr 
prooiems 


We  review1  nere  tne  BOC  aigonuim  in  soiling  me  svstem 

-lx  -  b 

ia)  Solve  .ix  -  b  using  tne  analog  opucai  processor  to  gei 

x„. 

ibl  With  a  dedicated  digital  electronics  processor,  read  x 
and  evaluate  tne  residue 

r,  -  b  -  -lx,  -  4\  -  -ix  =  -it  Ax„  l  (181 

ici  Normalize  a  to  use  ide  dvnanuc  range  of  the  svstem 
td)  Solve  opucaJIv  the  svstem 

■l  r  -  j  r  (191 

where 

t  «  it  Ax.,  i .  CO) 

and  s  is  the  radix  used  in  normalizing  r 
ici  Evaluate  eleciromcaJIv 

v  «  \  -  Ax  Cl) 

and 

nr .  m  b  -  -tv  C2) 

if)  If  r  is  small  enougft.  stop  Otherwise,  go  io  ici  and 
recvcie 

In  the  following  section  we  present  some  oi  me  preliminary 
results  from  computer  simutauon  studies  c:  me  BOC  in 
processing  adaptive  arrav  proDlems 

i.  SIMULATION  RESULTS 

Two  simulation  expenmenis  are  presented  :n  inis  section  in 
the  hrst  experiment  we  used  a  nve  element  arrav  and  assume 
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the  directions  of  the  lammcrs  were  Known  In  tne  second 
expenmen;  a  two  element  airav  is  used  ana  no  a  pnon 
information  is  assumed 

In  Figure  3  the  live  element  array  paitem  is  plotted  as  a 
luncuon  of  the  angle  v  Figure  3(a)  shows  tne  arrav  pattern 
before  adaptation  In  Figure  3(b)  the  pattern  alter  adaptation 
is  shown  for  a  jammer  at  45°  The  arrav  pattern  after  adapta¬ 
tion  nas  reformed  in  such  a  wav  that  it  nulls  the  lammer 
signal.  In  Figure  3(c)  four  jammers  are  considered  at  45°,  80°. 
120°.  and  150°.  The  atTav  pattern  is  again  reformed  to  null  all 
the  jammers  signal  reception 

In  Figure  4  the  BOC  was  used  to  solve  the  adaptation 
problem  assuming  no  a  prion  informauon  about  the  inter¬ 
ference  signals.  Figure  4<a)  shows  the  two-element  array  pat¬ 
tern  before  adaptauon.  In  Figure  4(b)-(d)  the  pattern  is 
plotted  for  a  single  jammer  placed  at  30°.  45°.  and  60°. 
respectively.  In  all  these  plots  the  arrav  adapted  to  cancel  the 
interference  signal  in  each  of  the  given  cases.  In  all  of  the 
preceding  results  the  tamrner  signals  are  considered  to  be  of 
the  same  strength  as  the  desired  signal,  and  the  convergence  of 
the  solution  obtained  in  less  than  five  iterations.  Also  the 
condiuon  numbered  of  the  /?,,  is  between  10°  and  oc. 


S.  CONCLUSIONS 

The  btmodal  opucai  computer  is  shown  in  these  preliminary 
results  to  present  a  powerful  mean  lor  solving  adapuve  phased 
array  problems.  We  are  considering  in  future  work  larger 
arrav  sues,  receiver  noise,  and  verv  strong  interference  signals 
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of 
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RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (CSI)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  CSI  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability / maintainability  and  compatibility. 


